FSU-CSIT-00-27 
JLAB-THY-00-34 



1 



Exact Chiral Symmetry for Domain Wall Fermions with Finite L s * 

Robert G. Edwards a and Urs M. Heller b 

a Jefferson Lab, 12000 Jefferson Avenue, MS 12H2, Newport News, VA 23606, USA 
b CSlT, Florida State University, Tallahassee, FL 32306-4120, USA 

We show how the standard domain wall action can be simply modified to allow arbitrarily exact chiral symmetry 
at finite fifth dimensional extent L s . We note that the method can be used for both quenched and dynamical 
calculations. We test the method using smooth and thermalized gauge field configurations. We also make 
comparisons of the performance (cost) of the domain wall operator for spectroscopy compared to other methods 
such as the overlap-Dirac operator and find both methods are comparable in cost. 



1. The problem with small eigenvalues 

Recently a great deal of theoretical progress has 
been made in the construction of lattice regular- 
izations of fermions with good chiral properties 
(HQ). For use in practical numerical simulations, 
though, approximations to these formulations are 
necessary. In the formulation using domain wall 
fermions (DWF) (l]||, the extent of an auxiliary 
fifth dimension has to be kept finite in numerical 
simulations while chiral symmetry holds strictly 
only in the limit of infinite fifth dimension. The 
violations of chiral symmetry are expected to be 
suppressed exponentially in the extent of the fifth 
dimension but in practice the coefficient in 
the exponent can be quite small |^,^| and the sup- 
pression correspondingly slow. 

In the case of overlap fermions there is no such 
problem in principle. However, there is a prob- 
lem of practicality: how to deal efficiently with 
e(H) — H/^/H 2 , where H is some auxiliary Her- 
mitian lattice Dirac operator for large negative 
mass, but free of doublers. Most commonly, the 
Hcrmitian Wilson-Dirac operator H w is used. 

Any good numerical approximation to e(H) 
must retain the property e(H) 2 = 1 which is cru- 
cial for the Ginsparg- Wilson relation 

{75, D ov {0)} = 2D ov (0) l5 D ov {0) (1) 

to hold and exact chiral symmetry for the overlap 



Dirac operator 

D ov {m) = ^[l + m+(l-m) l5 e(H)} (2) 

at in — to be preserved. A good approximation 
of e(x) is difficult to achieve for small x. 

Unfortunately, at least in quenched simula- 
tions with the Wilson gauge action, H W (M) = 
j 5 D w (-M), with D W {M) the Wilson-Dirac op- 
erator, has a non- vanishing density of zero eigen- 
values, p(0;M), for any gauge coupling j3 [[§. 
p(0; M) decreases quite rapidly with increasing 
/3, roughly as p(0; M)/cr 3 ^ ~ e" 6 " [§, but at 
currently used couplings a considerable number 
of modes with small eigenvalues exists. 

In all numerical methods to implement overlap 
fermions one can enforce accuracy of the approxi- 
mation to e(H w ) by projecting out the lowest few 
H w eigenvectors and adding their correct contri- 
bution exactly (see Ref. Q): with H w Vi — XiVi, 

n 

e(H„) = ^\v i )e()*){v i \+p£ ) App[e{H a )]pto ) 

i=l 

n 

Pi n) = l-$>^l- (3) 

i=l 

The number n of projected eigenvectors is cho- 
sen such that the approximation App[e(iJ„,)] to 
e(H w ) is sufficiently accurate in the subspace 
spanned by V j™^ . 
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2. Implications for domain wall fermions 

Integrating out the heavy fermions, one can 
show that the effective Dirac operator for the 
physical 4-d fermions from the domain wall 
fermion approach is pO|Jll[ | 

D t0 v(m) = - 1 + m + (1 - mh 5 E Ls/2 (H T ) (4) 

with £at(x) the polar decomposition approxima- 
tion to e(x) IJ) 



£ N (x) = 



and Fl 



(1 



\2N 



(1 



\2N 



(1 + x) 2N + (1 - x) 

tanh(2iVa;) 



2 A' 



H T {M)=H W {M)- 



1 5 H W (M) 



(5) 



(6) 



Zero modes of H W (M) are also zero modes of 
H T (M)\ 

Small eigenvalues of Ht induce chiral symme- 
try violations, unless L s is very large — tanh(j/) 
is within about 10~ 5 of 1 for y > 6.1. And the 
density of small eigenvalues of Ht is the same as 
the density of small eigenvalues of H w \ 

Manifestations of the induced chiral symme- 
try breaking are, for example, violations of the 
Gellmann-Oakes-Renner relation 



mXi = to 



;5>(a07r(O)) = <W) 



(7) 



finding a non-zero pion mass for massless quarks, 
and absence of exact zero modes in topologically 
non-trivial gauge field backgrounds. This latter 
is illustrated in Fig. [I] where we show the spec- 
tral flow of the hermitian domain wall operator 
H DW (Q;M) = l5 JD%(Q;M), with J the in- 
version operator of the fifth direction, the her- 
mitian Wilson operator H W (M) and the hermi- 
tian overlap H ov (0; M) = 75-D o «(0 : M) on a sin- 
gle instanton configuration. Right after a Wilson 
eigenvalue crosses zero the overlap operator has 
an exact zero mode while the lowest eigenvalue of 
the domain wall operator is clearly non-zero even 
for L s = 32 for M considerably after the crossing. 
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Figure 1. Eigenvalues of the hermitian domain 
wall, Wilson and overlap operators as function of 
the negative Wilson mass value M on a single 
instanton configuration. 



3. Exact chiral symmetry at finite L s 

With finite lattice spacing in the 5 th direction, 
0,5, the domain wall fermion (DWF) action is 



S DW = . 



(8) 



The kernel of the 5-d operator is, including the 
boundary conditions used to introduce an explicit 
quark mass to, 



£> (5) - 

DW " 

(D+ - iP_ 

-P + 



-P_ 
D + 



(9) 

(to - A)P+\ 





V mP- 





D+ 

-P+ 





-P_ 



where P± = (l±j 5 )/2 and D+ = a 5 D w (-M) + l. 
A(m) is the new ingredient, acting on the light 
fermions, chosen to insure chiral symmetry. We 
integrate out the L s — 1 extra fermion fields and 
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the pseudo-fermions to obtain |l] 
D tov (m;H T ) 



v- 1 



(5) 



1 + m + (1 - m)^£L B /2{a^H T ) 
term with A(m) . 



(10) 



Here En(x) is the polar decomposition approxi- 
mation to e{x) introduced in (||), and V is such 
that ('P _1 5 , )i = q is the light fermion field, 

T> = | P -^' & + P +^ + !. fe fOT 3 < L s (n) 

Jk \p_S Ls , k + P + S hk for j=L s . 

Ht, given in (^]), is the auxiliary Hamiltonian as- 
sociated with the transfer matrix T in the 5 th 
direction, 



T = 



1 - a 5 H T 



1 + a 5 H T 

With HtVi = \{Vi, Tvi — TiVi we set, 



(12) 



A(m) = {l-m) l5 {H w P_-l)Y j g l \v l ){v l \ , (13) 



9, = 2 



(rr L « - l) + (rr L ° + l) E (A,)] . (14) 



This "projects" low- lying eigenvalues of Ht from 
e L s /2{ a bHT) and gives the correct e(A;). Then 
Dtov{m] Ht) becomes D ov (m; Ht), 

D ov (m: H T ) = \ [1 + m + (1 - m) 75 e(ff T )] , (15) 
the overlap operator with auxiliary Hamiltonian 

With projection of the 5 lowest eigenvectors 
of Ht we now find the expected zero-modes 
of the massless projected domain wall opera- 
tor J (Dow(Q)V D^ w (0) for the same instanton 
configuration as in Fig. [l already for moderate 
L s , as can be seen in Fig. |. 

As already mentioned, the Gcllmann-Oakes- 
Renner relation (0) is a sensitive test of chiral 
symmetry. For overlap fermions 



m{b\{ lb b- v \m)f\b) = {b\b- v \m)\b) 



(16) 
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Figure 2. Eigenvalues of the domain wall operator 
with projection as function of the negative Wilson 
mass value M on a single instanton configuration. 
Also show is the Wilson operator spectral flow 
and the eigenvalues of 2Ht(M). 

holds for any chiral vector 75 \b) — ±\b), and aver- 
aging over several chiral Gaussian vectors \b) this 
becomes a stochastic estimate for the Gellmann- 
Oakes-Renner relation (Q). In ( |l6| ) D^(m) is 
the external fermion propagator with the con- 
tact term subtracted and multiplicatively normal- 
ized M, which avoids 0(a) effects pf, 



(l-m)D- 1 (m) = [D- v \m)-l] 



(17) 



D%{m) D^(1)V 



(5) 



- 1 
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We note that for standard DWF without pro- 
jection the subtraction and multiplicative nor- 
malization is done automatically, i.e., with q = 



(VJT)i where J tJ = <5 4 ,l s +i-j 
operator of the fifth direction, 

(qq)^{v- 1 [D^ w (^T 1 J'P} 11 

= T ^—[D^{m)-l] . 



is the inversion 



(18) 



Finite L s and no projection lead to strong vi- 
olations of the Gellmann-Oakes-Renner relation 
(0). With projection of a sufficient number of 
low eigenvectors the GMOR relation holds down 
to much lower quark masses, as seen in Fig. 0. 
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GMOR: j3 = 5.B5, B J x32, M=1.65, 50 cfgs 
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Figure 3. The ratio R = mx-x/{4>4>) for DWF 
without and with projection of the lowest 20 
eigenvectors. 



Figure 4. The cost of various implementations 
of chiral fcrmions as measured by the number of 
applications of D w on a 4-d vector. 
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4. Comparing Efficiencies 

Having different possibilities for chiral fermion 
formulations on the lattice we can compare their 
efficiencies in numerical simulations. For this we 
consider a spectroscopy calculation for a topologi- 
cally trivial gauge configuration. The cost is mea- 
sured in the number of times D w , common to all 
approaches, has to be applied on a 4-d vector. 
This is shown in Fig. ||. 

For a single quark mass preconditioned DWF, 
with projection to assure minimal induced chiral 
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spectroscopy calculation with overlap fermions 
one can use a multi-shift inverter to compute for 
all quark masses simultaneously, compensating 
for the factor 3-4 in increased cost. Which reg- 
ularization to prefer depends thus on the details 
of the numerical simulation attempted. If DWF 
are chosen then the projection method described 
here and in should definitely be used to en- 
sure good chiral properties at moderate L s . 
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